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ABSTRACT 


Receiver  operating  characteristics  for  phase-incoherent  detection  of  M 
multiple  observations  are  presented,  for  a  wide  range  of  values  of  M  and  signal- 
to-noise  ratio.  Comparison  with  a  processor  capable  of  coherently  summing  up 
all  M  signal  components  is  made.  The  additional  signal -to-noise  ratio  required 
for  the  phase-incoherent  combiner  is  found  to  be  greater  at  the  lower  false-alarm 
and  detection  probabilities. 
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INTRODUCTION 

In  many  cases,  received  signals  contain  multiple  (M-ary)  replicas  of  the 
signal  information  at  predictable  locations  in  the  time  and/or  frequency  domain. 
Knowledge  of  the  signal  replicas  can  be  obtained  by  analysis  in  a  noise-free, 
or  low  noise  environment,  by  placing  measurement  transducers  at  positions 
closer  than  is  practical  under  normal  operating  conditions,  or  by  utilizing 
artificial  means  to  force  higher  signal  levels  than  are  usually  available.  Once 
these  multiple  features  of  the  signal  are  known,  they  can  be  utilized  to  enhance 
the  detection  capability  of  a  system  operating  on  them.  However,  the  complexity 
of  any  system  that  attempts  to  detect  by  use  of  multiple  signal  characteristics 
will  be  increased. 

This  memorandum  examines  the  gain  available  when  phase -coherent  vs. 
phase-incoherent  detection  is  used,  as  a  function  of  the  number  M  of  signal 
characteristics  and  the  allowable  false  alarm  probability.  It  is  intended  to 
supply  the  information  required  to  make  the  tradeoff  between  increased  detection 
capability  and  system  complexity. 

The  analytical  results  in  this  memorandum  are  not  new;  for  example,  Helstrom 
[Ref.  1,  Ch.  Vl]  presents  all  the  equations  required.  However,  there  is  inade¬ 
quate  quantitative  evaluation  of  the  M-ary  phase-incoherent  processor  (Ref.  1 
contains  one  figure  on  page  181  for  a  false  alarm  probability  of  10  ).  Thus, 
comparison  with  an  ideal  processor  is  not  possible  for  a  significant  range  of 
signal -to-noise  ratios,  values  of  M,  and  false-alarm  probabilities.  The  present 
results  rectify  this  situation  and  make  the  comparison.  The  method  of  computation 
of  the  operating  characteristics  is  given  in  Refs.  2  and  3. 

PROBLEM  DEFINITION 

The  received  signal  is  composed  of  M  components,  each  deterministic 
(except  for  unknown  phase)  and  narrowband;  also  the  signal  components  are 
essentially  disjoint  in  time  and/or  frequency.  The  accompanying  noise  is 
zero-mean  stationary  Gaussian. 
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The  receiver  processor  to  be  considered  is  depicted  in  Fig.  1 .  Input  process 


Fig.  1.  Receiver  Processor 


x(t)  is  observed  over  time  interval  (a,b)  and  passed  through  a  parallel  set  of 
M  filters,  each  fairly  well  matched  in  time  duration  and  frequency  location  to 
one  of  the  M  components  of  the  received  signal.  The  filter  outputs  are 
envelope -detected  and  sampled  at  different  time  instants,  yielding  a  set  of 
outputs  (ijY} ,  •  The  method  of  combining  these  outputs  for  decisions  on 
signal  presence  will  be  discussed  shortly.  Notice  that  the  filters  [hjlrj] 
are  not  required  to  be  matched  filters,  but  rather  only  have  an  impulse  response 
duration  approximately  the  same  as  the  corresponding  signal  component,  and 
have  a  frequency  location  approximately  covering  the  signal  components' 
spectral  range.  This  allows  a  wide  variety  of  sub-optimal  receiver  realizations 
to  be  considered  and  compared.  Losses  in  performance  due  to  mismatched 
receiving  filters  or  improper  gain  settings  of  the  filters  are  easily  evaluated. 

If  just  one  signal  component  (M=l)  were  received,  and  phase-coherent 
reception  were  employed, instead  of  the  phase-incoherent  reception  technique  of 
Fig.  1  (i .  e. ,  no  envelope  detectors),  the  false  alarm  and  detection  probabilities 
are  derived  in  Appendix  A  for  an  arbitrary  filter  and  sampling  instant,  and  then 
specialized  to  various  optimum  filter  choices  for  different  observation  intervals 
and  noise  spectra.  On  the  other  hand,  for  phase-incoherent  reception  of  one 
signal  component,  the  false  alarm  and  detection  probabilities  are  derived  in 
Appendix  B,  under  conditions  similar  to  Appendix  A.  These  results  are  given 
for  the  sake  of  completeness  and  to  indicate  the  method  of  performance  comparison 
to  be  adopted  for  the  M-component  signal  case. 
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For  the  M-component  signal  case,  the  optimum  method  of  combining  outputs 
in  Fig.  1  (even  for  non-optimum,  but  given,  forms  of  filters 
is  derived  in  Appendix  C.  The  difficulty  of  implementing  and  analyzing  this 
processor  suggests  the  simpler  sub-optimum  processor  which  combines  the  outputs 
according  to 

.  0) 

j-' 

and  compares  z  with  a  threshold.  The  characteristic  function  of  decision 
variable  z  is  derived  in  Appendix  C,  and  the  false  alarm  and  detection  prob¬ 
abilities  given.  Specialization  of  the  general  filter  characteristics  to  optimum 
choices  is  also  presented  in  Appendix  C. 

If  it  were  possible  to  process  the  received  waveform  x(t)  in  Fig.  1  by  a 
phase-coherent  technique  (which  would  require  knowledge  of  the  phase  of  each 
signal  component),  the  performance,  in  terms  of  false  alarm  and  detection 
probabilities,  would  be  significantly  improved  for  large  M.  It  is  the  purpose 
of  this  memo  to  make  this  comparison  quantitative,  over  a  wide  range  of  false 
alarm  probabilities  and  number  of  components  M.  The  false  alarm  and  detection 
probabilities  for  coherent  reception  of  the  M-component  signal  are  presented  in 
Appendix  D. 


RESULTS 


The  characteristic  function  of  decision  variable  z  in  (1)  is  given  by  (C17), 
and  the  false-alarm  and  detection  probabilities  are  given  by  (C22)  and  (09) 
respectively.  The  performance  of  the  processor  in  Fig.  1  is  characterized  by 
the  single  parameter  dT,  which  is  a  measure  of  the  total  "  output  signal -to-noise 
(voltage)  ratio.  "  Specifically,  from  (03),  it  is  given  by 


<C  - 


^  ifrkjbi-^wr 

y B '  zff  dr,  dr,  Vi. (Tj  --cOK* (Tj  ’ 


(2) 


where  h^  (c)  is  the  complex  envelope  of  filter  hj  (c),  s_j(t)  is  the  complex 
envelope  of  the  j-th  component  of  the  received  signal,  and  RB(t)  is  the 
correlation  function  of  the  complex  envelope  of  the  received  noise.  The  filters 
in  (2)  are  not  assumed  to  be  matched  filters,  but  are  general,  subject  only  to 
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the  constraint  that  their  impulse  response  durations  and  bandwidths  be  approximately 
the  same  as  the  corresponding  signals.  The  detailed  assumptions  and  derivations 
leading  to  (2)  are  documented  in  Appendix  C. 

Plots  of  detection  probability  P„  (in  per  cent)  versus  dx,  for  false  alarm 
probability  PF  equal  to  10*  10*^  I CT*  and  10*  are  presented  in  Figs.  2-17,  for 
M=1  (1 )  10,  16,  32,  64,  128  ,  256  ,  512*,  respectively.  (There  is  a  change  of 
abscissa  scale  for  M-32).  As  M  increases,  the  performance  for  a  given  value  of 
dT  degrades.  For  example,  to  attain  PF  =  10  ,  PD=  .5,  a  value  for  dTof  6  is 
required  for  M=l,  whereas  a  value  for  c4  of  16. 6  is  required  for  M=512;  this 
is  an  increase  of  20  log(16.6/6)  =  8.  84dB  in  signal-to-noise  ratio.  Thus 
fractionalization  of  a  fixed  amount  of  signal  energy  degrades  the  processor 
performance. 

This  behavior  is  made  more  obvious  by  constructing  the  cross-plots  in  Figs. 

18-21  of  detection  probability  versus  dT ,  with  false  alarm  probability  PF 
held  fixed,  for  M=  1(1)6,8,10.  (The  curves  for  larger  values  of  M  can  be 
found,  if  desired,  from  Figs.  12-17). 

Figures  2-21  constitute  the  receiver  operating  characteristics  (ROC)  for  the 
processor  depicted  in  Fig.  1  followed  by  squaring  and  summation.  We  now  wish 
to  compare  this  performance  with  that  of  a  processor  capable  of  coherently 
summing  up  all  M  signal  components.  (We  do  not  necessarily  require  the 
coherent  processor  to  be  optimum;  this  aspect  of  the  problem  is  considered  in 
Appendix  D,  where  connections  with  the  optimum  filters  for  the  phase-incoherent 
processor  are  made).  The  false  alarm  and  detection  probabilities  of  the  coherent 
processor  are  given  in  (A13),  where  the  parameter  dc  is  given  by  (A12). 

The  method  of  comparison  of  the  processors  is  as  follows:  both  processors  are 
required  to  yield  the  same  PF  ,  Pft.  The  required  numerical  value  of  dc  for  the 
coherent  processor  is  then  smaller  than  the  value  of  d,-  required  for  the  phase- 
incoherent  processor  of  Fig.  1 .  The  increased  amount  of  signal-to-noise  ratio 
required  by  the  incoherent  processor  is  then  given  by 


*The  crosses  on  Fig.  17  will  be  explained  in  the  next  section. 
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20  Io^C^t/^c)  j  iBj  (3) 

since  both  dT  and  dc  are  measures  of  signal-to-noise  voltage  ratios.  (See  also 
Appendices  B  and  D  for  comparison  of  the  optimum  processors).  In  Figs.  22-24, 

(3)  is  plotted  versus  the  number  of  signal  components  M,  with  false-alarm 
probability  PF  as  a  parameter,  for  detection  probability  P*,  equal  to  .5,  ,9,  .99, 
respectively.  (There  is  an  ordinate  scale  change  in  Fig.  23.  Also  straight  lines 
have  been  drawn  between  the  integer  abscissa  values  for  ease  of  interpretation). 

It  is  seen  that  for  a  specified  detection  probability,  the  increase  in  required 
signal-to-noise  ratio  is  greater  for  the  lower  false-alarm  probabilities;  for 
example,  for  PD  =-5,  M=10,  the  increase  in  signal-to-noise  ratio  must  be 
5.33dB  at  PF=  102  but  only  2.69  dB  at  PF=  10  .  It  is  also  to  be  noted  that  the 
increase  in  required  signal-to-noise  ratio  for  the  phase-incoherent  processor  is 
greater  for  the  lower  detection  probabilities;  for  example,  for  PF=  10  ,  M  =  1  0, 
the  numbers  are  5.33dB  at  l^=  .5,  4.  15dB  at  PD=  .9,  and  3.55dB  at  PD=  .99. 

Thus  the  greatest  discrepancy  in  performance  between  the  coherent  and  incoherent 
processors  occurs  at  the  lower  detection  and  false-alarm  probabilities. 

In  Figs.  25-27,  similar  results  extending  up  to  M  =  512  are  plotted  on 
logarithmic  paper.  The  curves  appear  to  be  approaching  a  5  logM  dependence 
for  large  M.  For  large  M,  significantly  larger  signal-to-noise  ratios  are 
required;  for  example  1 3d B  more  is  required  at  Pp=  .5,  PF=  10,  M  =  512  for  the 
incoherent  versus  coherent  processor.  (Recall  that  the  ordinates  of  Figs.  22-27 
are  the  increased  signal-to-noise  ratios,  and  not  the  actual  values  of  signal-to- 
noise  ratio  required  for  a  specified  performance). 

GAUSSIAN  APPROXIMATION 

In  Appendix  E,  the  decision  variable  z  in  (1)  is  approximated  by  a  Gaussian 
random  variable  with  the  same  mean  and  variance.  The  detection  and  false-alarm 
probabilities  are  given  by  (E12).  For  large  M,  this  is  expected  to  be  a  good 
approximation,  since  (1)  is  the  sum  of  a  large  number  of  independent  random 
variables.  For  M  =  512,  the  Gaussian  approximation  is  plotted  as  crosses  on 
Fig.  17;  it  is  seen  to  be  a  very  good  one  for  PF=  10  ,  but  slightly  poorer  at 
PF=  10  .  Thus  the  Gaussian  approximation  can  be  used  for  other  values  of  M 
with  a  known  degree  of  accuracy  established  from  Fig.  17. 
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APPENDIX  A 

PHASE-COHERENT  DETECTION  (PCD) 

Consider  the  linear  filter  depicted  in  Fig.  A1 .  Input  process  x(t)  is 


Fig.  Al.  Linear  Filter 

observed  for  ,  and  passed  through  filter  h.  The  input  is  given  by 


Noise  nft>)  is  zero-mean  stationary  Gaussian,  with  correlation  1?^,  and  double¬ 
sided  power  spectrum  Signal  sfO  is  deterministic,  with  voltage  density 

spectrum  S(f).  We  have 


(Al) 


where 


y(t)=  ysttVynH:), 

i^(fc)=  Jo«h: 
yn(t)=  J^4-c 


(A2) 


(A3) 


Suppose  we  sample  output  «}#)  at  time  T  (which  is  arbitrary),  and  compare 
ytr)  with  a  threshold  L  for  purposes  of  signal  detection: 

2  *  y  (x)  X  L  . 

Random  variable  (RV)  i  is  Gaussian,  with 

0 

OR 


(A4) 


Etf  - 


and 


z}  =  II  d-x,  dxu  Ut-I,)Wt-xX  (V^)  s  cr\ 


(A5) 


(A6) 


a 
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The  probability  density  function  (PDF)  of  e  for  signal  present  is 
The  probability  of  detection  is  then 

J  ^  rM  = 

L 

where  normalized  threshold 

-A  *  L/T, 

and 

§(A '  il*  (?-)'*  «rK/2)- 

The  probability  of  false  alarm  is  obtained  by  setting  m*0  in  (A8): 

TF  =  5  i-X) . 

Now  define 

k 

cJ  -  J!l-  wtcw  Jt  -Pi Her  otfbyifc  _  ^ (T-w 

t_lr=  *****  oM  '  [fj<,r|,TiV,fr-T^fr-T^fc-^1 


(A7) 

(A8) 

(A9) 

(A10) 

(All) 


TV»en 


fF= 


35(A12) 


(A13) 


A  can  be  chosen  for  a  desired  TF  (Ref.  4,  Tables  26.5  and  26.6);  the  actual 
Tp  depends  additionally  on  dc .  Curves  of  performance  aregiven  in  Ref.  1,  Figs. 
IV.  1  and  IV.  2. 


Letting  u»^r  in  (A 7),  the  PDFs  for  signal  present  and  absent  are,  respectively 
fM  =  CM*  exp[-  i(u-d.T] ,  (A14) 

fo^  =(2^jVl  exp  [-  i  urji 

and  depend  only  on  the  parameter  d,  .  Thus  the  receiver  operating  characteristics 
(ROC)  depend  only  on  dc . 


Equations  (A1 2)  and  (A13)  are  the  most  general  situation;  filter  h  is 
arbitrary.  We  can  investigate  loss  of  performance  due  to  wrong  sampling  ig: 
T,  or  lack  of  knowledge  of  s  or  "R„,  or  any  realization  of  filter  h.  (d 
could  be  called  the  filter  output  power-SNR,  but  it  is  not  necessary). 
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Several  specializations  of  the  general  results  above  are  possible;  they  are 
considered  below. 


Case  1  Optimum  filter 

Choose  h  such  that  d  is  maximized.  It  is  the  solution  h#  of  the  integral 


equation 


f  4T^1(T-T)'RM(t-'r)  =  -sfO,  a<'ir<b. 


(A15) 


The  maximum  value  of  d  is  then  d  ,  where 


4*  -  J  dr  V)  (t-t)  sfr).  *  (A16) 

Case  2  Optimum  filter;  Infinite  Observation  Interval;  General  Noise  Spectrum 

As  a-*  —  «*»,  +°°  ,  the  integral  equation  (A15)  has  the  solution 

(A17) 


u  IL\_  ,  S*(f)  (r  i  2wf 


and  the  maximum  d  is  d  .where 
c  c 


(A18) 


—  °o  ^n\ 

Case  3  Optimum  filter;  Finite  Observation  Interval;  White  Noise 
For  white  noise,  but  finite  observation  interval, 

N(j| }  aN  ■p  ^(JouVle- sided 

•R.W-  Mj  Sto, 

and  the  integral  equation  (A15)  has  the  solution 

Vl3 Ct--c)  H<T<  1, 

with  a  maximum  d  of  d  : 


(A19) 


(A20i 
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K =  if  ’  <A2,) 

where  E  ,  is  the  received  signal  energy  in  (a,  b). 

Case  4  Optimum  Filter;  Infinite  Observation  Interval;  White  Noise 

The  maximum  value  of  d  is  d  .where 

c  c^ 

Ea®  E;  VT-T)  ®  T  i  K  =  >  (A22) 

and  E  is  the  total  received  signal  energy  over  all  time. 
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APPENDIX  B 


PHASE-INCOHERENT  DETECTION  (PID) 


The  processor  of  interest  here  is  depicted  in  Fig.  B1 .  Input  process  xft) 


xtt) 

Mr) 

»«)  r 

Eoy«.l  ope 

wfO 

De'tec'tov- 

Figure  B1 .  Linear  Filter  and  Envelope  Detector 


is  observed  for  fr  e 


and  passed  through  filter  h.  The  input  is  given  by 


C  nft  ' 

x  ft)  =  <  o* 


(Bl) 


Noise  n(b) 
Signal  sfr) 
Filter  hfc) 


is  zero-mean  stationary  Gaussian,  with  narrow-band  spectrum  G  . 
is  deterministic  except  for  unknown  phase, and  narrow-band, 
is  narrow-band. 


Let* 

5  It)  «  “Re  { 1  ft)  0  2-rrii+  i  &)}  > 

n  ft)  =  'Re  £«  ft)  exp(i  ) 

\\  tc)  =  Ke  £h  he)  exp  (i  2^,"^ , 


RV  0  is  uniformly  distributed  over  2-tt .  The  same  center  frequency  f0  is  no 
restriction,  as  the  complex  envelopes  (underlined  quantities  in  (B2))  can  absorb 
differences  in  center  frequencies.  The  complex  envelope  of  is 


Complex  notation  is  used  freely  in  this  Appendix;  see  Ref.  1,  Ch.2, for  example. 


38 


Tech.  Memo.  TC-1 79-71 


y_  Hr)  s  J  -  £  dr  K  (t-r)[ s  (x->  e  &  +  «  ^)]  ^  Jj.,  ^  if  sign*!  p  r*s«ft .  (B3 ) 


The  output  of  Fig.  B1  is  given  by 


ww-IsM 

(B4) 

Suppose 

threshold 

we  sample  wtt)  at  time  T  (arbitrary),  and  compare  with  a 

L  for  purposes  of  signal  detection: 

Zrs  wCt)  >  L. 

(B5) 

Then  using 

(B4)  and  (B3^ 

4  .  |  +■  J».w| 

(B6) 

Let 

JfntT)  =  x  -h  i  j  . 

(B7) 

x  and  y 

are  real  Gaussian  RVs,  with 

(B8) 

And  since  (Ref.  1,  Ch.  2) 

*  0  - 

_  h  .  ^  (B9) 

|vjnCr^|a  =  4-  JJdxi  Jt,  )i  tr-x^  (x-T-^l^  [x-rx^  =  x  ty  j 

there  follows 

*3  = 
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Then 

yUT)  =  i|>wl  =  <r\ 

(BID 

Let 

=  "R+il  "fw 

(B12) 

Then  using 

(B6),  (B12),  and  (B7), 

(B13) 

|K^.\X+-x+\u1  =  |AViB|  = 

where 

A  h  'R+*) 

Then  the  PDF  of  real  RVs  A  and  B,  for  a  given  value  of  ©  ,  is 

The  probability  of  detection  for  this  given  value  of  6  is 


(©)  -  >  L  |  9) 


/S?>L 


(B14) 


(B15) 


(B16) 


If  we  substitute  (B15)  in  (B1 6)  and  let  A^cos^,  B  =  rsin4>, 
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•I*. .  ■  q(^2!.  i). 

7<r 

which  is  independent  of  B  .  (This  result  is  related  to  that  in  Helstrom,  Ref.  1, 
p.  154,  under  (3,16),  but  is  more  general  here  because  h  is  arbitrary). 


(B17) 


Letting 


JV  =  U<r, 


4 

we  have 


i  dr  i  (t-t)  s  (x)| 


(B18) 


[2  k^r)h\r~^9(rrz^j,/2  ’ 


V^x>^\ 

?F  -  Q 


(B19) 


can  be  chosen  for  a  desired  Pp  ;  the  actual  7?,  depends  additionally  on 
dz  .  Curves  of  performance  are  given  in  Ref.  1,  Figs.  V.2  and  V.3. 

If  we  let  u  =  2/<r  in  (B13),  it  follows  from  (B 1 4)  -  (B18)  that  the  PDFs 
for  signal  present  and  absent  are,  respectively, 


■p^u)  =  U  U>  0, 

p0(»*)=  u  o<p[-^t/);  U>0, 


(B20) 


and  depend  only  on  the  parameter  The  quantity  dE  has  the  following 
interpretation: 

2 


.  HjfcWL 

2.  rt-1 
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instantaneous  -PtBaer  s^nal  Pufedt.  p^uier  qt  -time  T 
"ihstatrttfweaus  filter  mo» sc  output  power  ut  fiwie  T 


(B21) 


The  denominator  of  (B21)  is  the  ensemble  average  filter  output  noise  power 
at  time  T.  Steady  state  need  not  have  been  reached  for  either  the  signal  or 
noise.  </2  could  be  called  the  instantaneous  filter  output  power  -  SNR  if 
desired;  this  is  not  necessary. 


(As  a  special  case  of  this  general  relation  for  dr  ,  consider  a  sine  wave  of 
amplitude  A  in  the  presence  of  noise,  and  an  observation  interval  long 
enough  for  the  filter  output  to  reach  steady  state.  Then 


£&)  *  A, 


(B22) 


and 


1  (.t)  =  j^tC  =  2  2  T*  *  Htr  eutprft  biwe  ujave  (B23) 


The  filter  output  noise  power  is 

yptj  =  <r". 


Than 


JLp* 

~2~  =  — ■  t  ■  st  -pilfaer  output  povJtr-SNR. 


(824) 


(B25) 


Several  specializations  of  the  general  results  (B18)  and  (B19)  now  follow: 
Case  1  Optimum  Filter 


Choose  h  such  that  <4  is  maximized, 
integral  equation 


It  is  the  solution  h»  of  the 


£ 


dr  ht  (t-x)  =  §  H),  a <  -fc<  b. 


(B26) 


The  maximum  value  of  Ar  is  then  ,  where 

I 


2|  dr  H,(t-t)5^, 


(B27) 
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(which  is  always  real  and  non-negative).  The  reason  this  result  for  <JI( 
differs  from  Helstrom,  Ref.  1,  p.  151,  (3.6)  by  a  factor  of  2  is  due  to  a 
different  definition  of  "Ra  .  Looking  at  Helstrom,  Ch,  II,  sect.  5,  we  find 
his  0(X)  =  a'Katx)  .  Therefore,  p.  147,  (2.16)  gives  QW  =  2)JLT-t), 
and  his  (3.6),  p.  151  becomes  identical  to  our  . 

Case  2  Optimum  Filter;  Infinite  Observation  Interval;  General  Noise  Spectrum 


As  a-*  —  °o,b~>  +  00,  the  integral  equation  (B26)  has  the  solution 

, ,  ( *  i  *0f)  exf  (-« Zirf'rt)  \  (B28) 

Dxlf;  = - sntr - 


and  the  maximum  ^  is  ,  where 

lS(f)P 


Grtf) 


(B29) 


Case  3  Optimum  Filter;  Finite  Observation  Interval;  White  Noise 

For  white  noise,  but  finite  observation  interval, 

Grw(4)  =  Nj  >  all  -f, 

,  all  I, 

'R,,  4  &  (v), 

and  the  integral  equation  (B26)  has  the  solution 


(B30) 


(B31) 


with  a  maximum  of  <1^  : 


(B32) 
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where  is  the  received  signal  energy  in  (a,b).  This  relation  for  <1x3  / 

combined  with  p((u)  in  (B20),  agrees  with  Davies,  Ref.  5  (remember  that 

N0  =  2^). 

Case  4  Optimum  Filter;  Infinite  Observation  Interval;  White  Noise 
The  maximum  value  of  d*  is  <1x4  /  where 

(t~t)  =  afl  T>  (B33) 

=  e/No 

and  E  is  the  total  received  signal  energy  over  all  time. 


Notice  from  (A21),  (A22),  (B32),  and  (B33)  that 

dCj  =  >  dc^  -  dl4  (B34) 

for  the  same  signal  waveform  and  noise  spectrum.  These  are  special  cases  of 
the  more  general  relation 

V  Ji. ,  o»5> 

which  holds  for  narrow-band  signals  and  noise.  To  see  this,  let 
I  li)  -  fa  £  S  Hr)  exp(i  , 

K,(r)  =  'Ke  (B36) 

"K», W  =  2^e^a^)exp0^iT^  Casing  (33^)), 


in  Case  1  of  PCD,  Appendix  A.  Using  the  high-frequency  behavior  of  <xp(i3»r{,i\ 
the  integral  equation  in  PCD,  Case  1,  becomes 


-i  in-^T 
e 


•RQ(-t_T)= 


(B37) 
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Defining 


„  a  i2irtT 

k(T-t)=fiMe  , 


(B38) 


this  last  integral  equation  becomes 


fVriT-x)v-M  =  a<t<  i>. 

Then  using  the  narrow-band  behavior  of  s  and  h,  again,  ic* 

4*  *  dr  i  })(T-T)e'^r  ihr)+-  ct>*ju^ 


becomes 


(B39) 


=  t^dx  ]\  (T-t)  -s  fx)  +■  <-0«*J>Wx  covj^a+e 


(B40) 


*  2] 


this  is  real  (remember  k  is  the  solution  of  integral  equation  (B39)).  But  now 
(B39)  and  (B40)  are  identical  to  Case  1  of  PCD,  upon  identifying 


k  (T-T)  =  )],  (t_t). 


(B41) 


or  equivalently 

J  (T-T)  4-  k,  (t-  x)  e 


or 


k  fr)  -  1  ii,  (r)e 


-1  2irk>T~ 


(B42) 


(B43) 
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Therefore 


dL  =  aT  (B44) 

ci  « 

for  a  given  signal  waveform  and  noise  spectrum.  Also,  the  optimum  filters  are 
identical  except  for  an  irrelevant  phase  in  PID.  (Scale  factors  are  unimportant). 
See  also  Helstrom,  Ref.  1,  p.  154,  for  a  similar  statement. 

However,  it  is  n&t  necessary  that  dc  equal  <lr  ,  as  can  be  seen  by 
sampling  y^ft)  in  PCD  at  a  zero  of  the  high-frequency  waveform;  then 
dc  =  0  while  dt  ^  0  .  In  fact,  d  can  be  negative,  whereas  for 

all  sampling  instants. 

Although  dI(  for  a  given  signal  waveform  and  noise  spectrum, 

performance  is  poorer  for  PID  than  for  PCD,  because  the  ROC  for  PID 
is  governed  by  the  Q-function,  whereas  the  ROC  for  PCD  is  governed  by 
the  $  -function. 

Suppose  instead,  we  want  PCD  and  PID  to  yield  the  same  performance, 
i.e-,  identical  Tp  .  Then  reading  the  corresponding  ROCs,  suppose  we  need 
values 

<U  ,=  V< 

V-vrt>vO  f»  Pi1*-  <B45) 

Then  the  increased  SNR  required  by  PID,  in  comparison  to  PCD,  is,  in  dB, 

20  loj,.(Vx/Vc).  (B46) 

(The  factor  20  is  due  to  the  fact  that  d  is  proportional  to  signal  voltage,  not 
power,  in  both  PCD  and  PID). 
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APPENDIX  C 


PHASE-INCOHERENT  DETECTION-  M  COMPONENTS  (PID-M) 


Input  process  xft)  Is  observed  for  t£  (a,b),  where 

C  ntfr) 

xft)  =  - S  OB  {  * 

where  ifct)  and  n(t)  are  narrowband  processes.  Input  signal  stt) 
of  M  components: 

M 

Sk)=  ^  Sv  It), 

KM 


(Cl) 


is  composed 


where 


s K(t)  -  ft)  exp(i2irit +i^ ;  sK  U)  =  Ak  Fk  (t). 


(C2) 


Ak  is  the (real) amplitude  of  the 
equal  amplitude.  RVs  {©kY^ 
over  a  2ir  interval . 


signal  component;  ^ft)],  all  have 
are  independent  and  uniformly  distributed 


Now  consider  the  following  processor: 


Vl,  (TO 


1 

ETivdoj?* 

w./t)  j 

5«mpie 

i 

J 

L>«st :«c  t»** 

d  r, 

x(t> 


i _ 

envelope 

f t)^ 

5«mpk 

«tT„ 

Figure  C.l.  Receiver  Processor  for  M  Components 
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The  filters  {V'jte)!  are  arbitrary  for  the  moment.  The  complex  envelope  of 

$!*■>  is 


M  i©k 
.e 

k»i 


+-  t  jNr  Lfc-r)  Jlfr). 

ft  J 


(C3) 


We  sample  Wjft)  at  time  Tj  .  The  signal  component  of  is,  forgiven 

{A}  - 


1  ift< 

k-\ 

We  shal  I  assume  that 


h  rb  . 

fft  ^  Kj  ^  s  Ak  ^  ^Tj  Fx  W  *  ©  ^  0  *  *■ 

This  means  that  the  j-th  branch  of  the  processor  does  not  respond  at  time  "Tj 
to  the  K-"A  input  signal  component  if  j  £  K  .  This  occurs,  for  example,  if 
the  signal  components  are  disjoint  with  each  other  in  time  and/or  frequency 
and  if  each  filter  is  cpproximately  matched  to  its  corresponding  signal;  this 
is  a  realistic  situation.  The  signal  component  of  ij_j  (ij)  is  then 

e'6j  SjW  -  Cj. 

The  noise  component  of  is 

nj  ‘ 

This  is  a  complex  Gaussian  RV  with  "n.  =  0  ,  since 
Also 

h.  \  ~  0,  Smce  ja.(*t,)  ia^k)  *  o. 

J 


(C4) 


(C5) 


(C6) 


(C7) 
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And  «jD*  - 

•  i  ^  <hi  *;  kj  h  k*  C-Tj  -^)Ka fc.-^)Sp , 

A 

where  we  have  used  assumption  (C5)  above.  Therefore,  the  noise  components 
in  are  all  independent. 

The  PDF  of 

^(Ti)  =  lco+nil 

is  available  from  (B6)-(B20)  as 

?>j>  =  ex|>[-  2 -^4 - i ,  V  °> 

j  «* 

Foh)  =  **f  [-  i  ^r]  ,  V0' 


(C8) 


(C9) 


(CIO) 


where 


r?*i  -  -J-  JJ  <n:,<te  ki(Ti-x-)k*('ii-x.)^(T;--ci)) 

^  A 

j  t  A  *  i£iL  ,  iCJlkJ(ij--^5itr)| 

XM/  r-\  w  . 


(Cll) 


Or  jf  <Mx*  kjCTJ-T^hj  (Tj-xO'Rfllrr^  a 


Due  to  the  independence  of  the  noise  components  in  ,  the  likelihood 

ratio  based  on  observations  is 

i  -  TT 

j=t  J 

-  ^p[“i  ^r]  Xs  (^xCi)lr  ), 

J*  ^ 


(Cl  2) 
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where 


l  1—1 


Cdr  kj(T;-r)^lx-)[ 


J  =  > 


i  5  <*,  Jta  h..  (lj  -T,)V)*  (-Tj  M*) 


The  optimum  processor  is  therefore  based  upon  a  comparison  of  the  RV 

M 


zi* 


with  a  threshold.  But  this  processor  requires  knowledge  of  absolute  signal 
levels  for  its  realization. 

We  instead  consider  the  sub-optimum  processor  (^uickJ  hjj  fcl4)) 

M 

z  =  Xwj 

j- j 

and  inquire  into  its  exact  performance.  From  the  PDF  in  (CIO),  the 

characteristic  function  (CF)  of  v£  for  signal  present  is 


(C13) 


(C14) 


(Cl  5) 


(Cl  6) 


At  this  point  we  will  assume  that  all  the  T?  are  equal.  Recalling  the 
definition  of  O'j'  in  (Cll),  this  is  tantamount  to  assuming  the  noise  has 
substantially  the  same  spectral  level  in  the  region  of  each  of  the  filters  {hj^  . 
(The  analysis  could  be  extended  to  unequal  ;  however,  the  computational 

effort  would  be  significantly  greater. )  Denoting  the  common  value  by  <T^Z  , 
the  CF  of  i  in  (Cl  5)  is 


(Cl  7) 


which  depends  on  only  through  if*  as  defined  in  (C13).  The  PDF  of 

2  is  (Ref.  1,  p.  174,  (2. 10)  and  (2.12)): 
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_! 


-■  (  jfr/)  “l*[-  i  ^  -  i  rM-,  ( s  a  0 


(Cl  8) 


Comparison  of  Z  with  a  threshold  L  yields  (Ref.  1,  p.  176) 
P  =  Prob(?>L)  = 


(Cl  9) 


where 


.  .  J£_ 

~  <rz  * 


(C20) 


and 


Also 


Qm  («,  (!»)  =  x  fc)  exp[*  Tx*~  TM-i  (<x). 
?F  =  QjM-')  = 

j=o  w  ’ 


(C21) 


(C22) 


The  performance  of  this  suboptimum  processor  depends  on  the  individual  signal 
strengths  only  through  the  quantity  (C13): 

4  =  5.  <C23> 

j-1 

Thus  it  does  not  depend  on  how  the  signal  energy  is  fractionalized  into  its 
individual  components;  ^  alone  counts.  This  holds  even  though  the  filters 
tw  in  Fig.  Cl  are  not  matched  to  the  received  signdl  components. 

The  performance  of  the  PID-M  processor  in  Fig.  Cl  followed  by  a  summation 
of  squared  outputs  can  be  optimized  by  maximizing  the  quantity  dT  .  This 
in  turn  can  be  maximized  by  choosing  filters  $.hj(r)i  such  that  [d^ij  )}  are 
maximized.  Let  be  proportional  to  the  optimum  complex  envelope  for 
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the  filter  which  maximizes  djCk)  at  time  when  is  received. 

Then  from  PID,  Appendix  B, 


(C24) 


Notice  that  these  filters  do  not  require  knowledge  of  the  relative  signal  strengths 
of  each  component.  The  maximum  for  each  component  is  d^(K)  /  where 


^  (k)  =  2]  dr  A  ktK  (T*-r)  sk  It)  =  A*  2  ^  <hr  k*  J^Crr),  (C25) 


and  we  have  used  (Cl  1),  (C24),  and  (C2).  The  entry  to  Tp  in  (C19)  is  d-^  ? 
where 

d-N  :S.dx((k).  (C26) 

*  K— t 


We  do  not  have  to  assume  infinite  observation  interval  or  white  noise  to 
reach  these  conclusions.  The  major  assumptions  are  the  disjoint  signals  (under 
(C5))  and  the  equal  variances  (under  (Cl 6)).  , 
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APPENDIX  D 

PHASE-COHERENT  DETECTION-M  COMPONENTS  (PCD-M) 

We  limit  consideration  to  the  optimum  receiving  filter  here;  Appendix  A 
covers  the  more  general  case.  The  phase-coherent  detector  collects  all  the 
signal  energy  of  the  M  components  and  yields 

pp- fK-A),Pp- <D,> 


where  (from  Appendix  A), 

dc*  -  s(t), 

1  A 

and  Ji,  is  the  solution  of  the  integral  equation 

f  dTh,(r--r)K^-T)  =  *(t),  a  <t<^  i>. 

s(t)  is  the  total  received  signal,  and  T  is  the  sampling  instant  (which  is 
arbitrary). 

Now,  from  (C2), 

S  Hr)  -  «.xj>(i2v£t  +  i9K)}. 

k- 1  '■ 

Let 

A 

j  K,KCT-T)'Rlllt-X)  =  TOT  IS 

Then  the  filter  ht  in  (D3)  is  given  by 

N\_ 

V»,(r-r)  -  KkCt-x), 


(D2) 


(D3) 


(D4) 


(D5) 


(D6) 
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Also  then  M  M 

<  =  f  k  dr  K  It-t)  5  (r) .  j  Jr  2  k.  tot)  2  s, 

L  (D7) 

M  *  A 

=  5.  2  *  k(r-to  ij  to 

K--  I  j  =  »  * 

Now  let  us  relate  the  of  this  section  to  the  £.h,K/T)}  in  (C24). 

To  this  end,  (using  narrowband  behavior  of  the  optimum  filter^  let 


l1K  (r)  =  Ke[h,KW  **)>('  ^r)j. 


(D8) 


Then  (D5)  becomes,  using  the  narrow-band  behavior  of  the  various  functions, 

A  “  ,  ie* 


[  dx  '  =  AKFK(t)e  , 


a<~t<  k. 


(D9) 


Comparing  this  with  (C24),  we  see  that 

•i"k,K^T'T')e-'^T  =  «<<r<  V. 


(DIO) 


Thus,  the  K-'tk  phase-coherent  filter  hiK  requires  knowledge  and  use  of 
whereas  the  phase-incoherent  filter  does  not.  Then  there  follows 

J  dT  V|,  (t- t)  S:  hr)  -  J  dr  i  jk  *  (r-r) e  lk  e  J 

a  J  « 

■^Ajtjtoe' Aj  Fj*^^  eihrf'r} 

-  ^  [  dr  k (K  (t- r) J  +-  c&wt,l*x  conjugate  (Dll) 


.  b  •  0  I Q. 

=  j  dt  AK  e. 1  K  h  ,K  tTK_r)Aj  F,(r)t  J  +-  coylu  lonjujait 


(D12) 
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=  An'  f  dr  i  *  (V^  Fk  hr)  £jK  +  <•*»?!«  ct”l)^e 


(D13) 


-  ijK  • 


(D14) 


The  transition  to  (Dll)  utilizes  the  high-frequency  character  of 

(D12)  utilizes  (DIO),  (D13)  utilizes  (C5),  and  (D14)  utilizes  (C25).  Then  from 

(D7)  and  (C26), 


M 


=  dr,  . 


(D15) 


Therefore 


(D16) 


That  is,  for  a  given  signal  waveform  and  noise  statistics,  dc,=  dTi  under 
the  assumptions  above.  However,  the  ROC  for  PID-M  is  poorer  than  the  ROC 
for  PCD-M,  because  the  former  is  governed  by  the  QM~  function,  whereas 
the  latter  is  governed  by  the  function. 


If  PID-M  and  PCD-M  are  to  yield  the  same  } 

dt,  -  Vc  for  PCD-  M, 

dT)=  Vr  (^Vc)  for  Plb-M. 


suppose  we  need  values 


(D17) 


Then  since  each  "d"  is  proportional  to  signal  voltage  (not  power),  the  increased 
SNR  required  by  the  PID-M  is,  in  dB, 


20  loj*  (  vx /vcy 


(D18) 
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APPENDIX  E 


GAUSSIAN  APPROXIMATION 

The  CF  of  decision  variable  2  in  (05)  is  given  by  (07).  The  CF  of 
the  normalized  variable  U  =  is  then 


r'2f)  . 


(ED 


For  large  M,  only  small  values  of  £  lead  to  significant  values  of  the  CF 
(El).  Then  since 


-i2f  s  txp^-i  2f  ^  4  order  f  > 


(E2) 


and 


5  *XF  l^(l+,^)dTj  ^  0,rcler 

(El )  becomes 

(^t  + 2M)  ~  order  T, 

The  PDF  corresponding  to  CF  (E4)  is 

t  r,(x- 

^t2  <«+*)  j  ’ 

which  is  a  Gaussian  RV  with  mean  and  variance 

(E5)  is  consistent  with  considering  the  RV 

M 

a  -  Z*; 

K  = » 


(E3) 


(E4) 


(E5) 


(E6) 
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as  a  Gaussian  RV.  The  statistics  of  are  available  upon  consideration 

of  (Cl 5)  and  (CIO);  namely 


X~J  -  2+  ,  \far txsi  =  -t[l+di(‘)]. 


(E7) 


There  follows  immediately  from  (E6),  (E7),  and  (C13), 
a  =  2fA  +  dr  ,  Vav-^u]  =  4(m+«C). 


(E8) 


The  probability  of  detection  follows  from  (E5)  as 

V  =TWl(»>T)  = 


where  T  is  a  threshold,  and  $  is  given  by  (A10).  The  false-alarm  probability 
is  obtained  by  setting  dr*  o: 


(E10) 


If  we  define  a  new  threshold  -A.  according  to 


(E9)  and  (E10)  become 


pF  = 


l\|d^+K  /• 


(Ell) 


(El  2) 


Equation  (El 2)  is  plotted  as  crosses  on  Fig.  17.  The  new  threshold  .A.  can  be 
chosen  for  specified  Pp  .  . 
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